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Abstract 

We study in detail the quasinormal modes of linear gravitational perturbations of plane- 
symmetric anti-de Sitter black holes. The wave equations are obtained by means of the Newman- 
Penrose formalism and the Chandrasekhar transformation theory. We show that oscillatory modes 



. decay exponentially with time such that these black holes are stable against gravitational pertur- 

cn : 

bations. Our numerical results show that in the large (small) black hole regime the frequencies 

m 

J> ■ of the ordinary quasinormal modes are proportional to the horizon radius r+ (wave number k). 

^3 . The frequency of the purely damped mode is very close to the algebraically special frequency in 

' the small horizon limit, and goes as ik'^/3r^ in the opposite limit. This result is confirmed by an 



analytical method based on the power series expansion of the frequency in terms of the horizon 



O . radius. The same procedure applied to the Schwarzschild anti-de Sitter spacetime proves that the 

$-H . purely damped frequency goes as i{l — l)(/-|-2)/3r+, where / is the quantum number characterizing 

the angular distribution. Finally, we study the limit of high overtones and find that the frequencies 



^ ■ become evenly spaced in this regime. The spacing of the frequency per unit horizon radius seems 

to be a universal quantity, in the sense that it is independent of the wave number, perturbation 
parity and black hole size. 
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I. INTRODUCTION 



The dynamical response of black holes to any external disturbance is dominated, at inter- 
mediate times, by a discrete set of damped oscillations called quasinormal modes (QNMs). 
This result has been verified not only at linearized level, but also in fully numerical simula- 
tions for stellar gravitational collapse jl | and collision of two black holes P| . The possibility 
of observation in the forthcoming years of the radiation due to the QNMs opens a new per- 
spective for the gravitational wave astronomy. From these observations it will be possible 
to obtain direct evidence of black holes' existence, as well as to estimate their parameters 
(mass, electric charge, and angular momentum) j^. 

In addition to the astrophysical applications, an interpretation appeared some years ago 



of the QNMs in terms of the AdS/CFT conjecture j^. According to this conjecture, there 
is a correspondence relation (or duality) between string theory in an anti-de Sitter (AdS) 
spacetime and a conformal field theory (CFT) on the boundary of this space PjO]. It implies 
that a large black hole in the AdS spacetime corresponds to an approximately thermal state 
in the CFT. Perturbing the black hole is equivalent to perturbing the thermal state, and the 
perturbation damping time gives the time scale for the return to the thermal equilibrium. 
Therefore, by computing quasinormal (QN) frequencies in an AdS space it is possible to 
obtain the thermalization time scale in the strongly coupled CFT. 

The early works studying perturbations of AdS black holes considered the evolution of 



a conformally invariant scalar field in such background spacetimes After the advent 

of the AdS/CFT conjecture, QNMs of spacetimes with a negative cosmolo gica l constant, 



such as Kottler (also callec 
Nordstrom-AdS ll 
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Schwarzschild-AdS) 
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Reissner- 



were computed for different 



no 



and Kerr-AdS 

fields, dimensions, and boundary conditions. Within this program, other kinds of black 
solutions were also analyzed. For instance, in the Banados-Teitelboim-Zanelli solution 
analytical expressions for the QN frequencies presented exact agreement with the poles of the 
retarded correlation function in the dual CFT 0, • This result provided an important 
quantitative test of the AdS/CFT correspondence. As a complement to the aforementioned 
studies, analyses of the time evo 



have been performed 
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ution of some fields in asymptotically AdS backgrounds 
IsO^ . More recently, the conjecture relating the highly 



damped QNMs to black hole area quantization has motivated a search for QN frequencies 
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with large imaginary parts also in asymptotically AdS spacetimes |3ll. l32l l33l. l34l l35l l36|. 

The studies of perturbations in AdS spacetimes were also extended to black holes with 
nonspherical topology. In Ref. jsl] the evolution of a conformal scalar field in the family 
of topological backgrounds was investigated, which includes the plane-symmetric black hole 
as one of its members, while in Ref. j29| the case of a nonminimally coupled scalar field in 
the same kind of spacetimes was considered. Lately, the lowest QN frequencies associated 
with the scalar, electromagnetic, and grav itational perturbations of the plane-symmetric 
black hole background were computed [37|. The study of the scalar and electromagnetic 
fields was relatively complete. However, important questions were left without answers in 
connection with the gravitational perturbations. For instance, what are the wave equations 
governing general metric perturbations? What is the physical interpretation of the modes 
with different wave numbers? What are the QN frequencies in the intermediate and small 
black hole regimes? What is the relation between these frequencies and the wave number? 
Why is the pure imaginary frequency not proportional to the horizon radius in the large 
black hole limit as it happens for ordinary modes? Is there some relation between the latter 
and the algebraically special frequency? And what is the behavior of the QNM spectrum in 
the regime of highly damped overtones? 

The objective of the present work is answering most of the previous questions as well 
as other ones that appeared during this work. Some of the conclusions drawn here are 
extended to the Schwarzschild-AdS black hole. The plan of this paper is the following. In 
Sec. m we write the equations defining the background spacetime in a Newman-Penrose 
(NP) null tetrad frame j^Sl- In Sec. IIIH we obtain Schrodingerlike equations governing the 
axial and polar perturbations. In Sec. IIVI we define the QNM boundary conditions to be 
used and prove the black holes stability against perturbations with a non-negative potential. 
A brief description of the Horowitz-Hubeny method (Sec. |V|) follows, as well as subsequent 
exposition of the QN frequency values for different horizon radii and wave numbers (Sec. 
IVIjl . Then, we study the large black hole limit of the axial QNMs (Sec. I Villi and the non- 
QN character of the algebraically special frequency (Sec. IVIIIIl . Finally, we conclude by 
summarizing and commenting the main results (Sec. IIX|) . 

Throughout this work, we use the geometric system of units in which the speed of light 
and the Newton's gravitational constant are set to unity. The background spacetime, more- 
over, has a natural lengthscale given by the AdS radius R = a/3/|Ac|, where A^ is the 
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negative cosmological constant. Here we choose the length unit so that R = 1 and measure 
everything in terms of the AdS radius. In relation to notation, the asterisk (*) denotes 
complex conjugation and the definition of the NP quantities follows Ref. js^. 

II. THE BLACK HOLE SPACETIME 

With the aim of investigating the evolution of linear gravitational perturbations in the 
background of a plane-symmetric AdS black hole [4^, we summarize here the basic properties 
of such a spacetime. The metric may be written in the form 

ds^ = f{r)de - f{r)-'dr^ - r\dip^ + dz^), (1) 

where 

/(r) = - 4M/r. 

There is an event horizon at r = r_|_ = (4M)^/^, and an essential singularity at r = 0. 
Depending on the ranges of coordinates ip and z, the metric may represent a cylindrical 

{0 < (p < 271, —00 < z < +00), toroidal {0 < ip, z < 27r) or planar (—00 < ip,z < +00) black 



hole spacetime. For the torus, M is related to the ADM system mass |41[; for the cylinder, 
to the mass per length unit of z line 0, Q|! plane, to the mass per area unit of 

the Lf — z plane |4J| . 

According to the Petrov classification, the above metric (^J is algebraically type D. It is 
then convenient to define the NP quantities in terms of a Kinnersley-like null frame |j45| . In 
this case, the two real vectors of such a frame, D and A, are chosen to lie, respectively, on 
the ingoing and outgoing radial null geodesies and are defined by 

D = Fd, = ji9t + fdr) , (2) 

A = n^d, = ^{dt-fdr), (3) 

which are the double principal null directions of the Weyl tensor. In order to complete the 
frame basis we have the vector 

6 = m^d, = ^{d^ + id,), (4) 

and its complex conjugate 5* = m*^d^. Both of these vectors, and m*^, are orthogonal 
to and ?7.^, and satisfy the normalization condition m^m*^ = — 1. With the chosen tetrad 



basis, the only nonvanishing spin coefficients are 



P 



1/r, 



//2r, 



7 = (r^ + 2M)/2r^ 



(5) 



and the only nonvanishing Weyl and Ricci scalars are 



^2 = -2M/r^ 



A = 1/2. 



(6) 



III. THE WAVE EQUATIONS 

Here the equations for the linear gravitational perturbations are obtained by using the NP 
formalism, which allows us to investigate general properties of such perturbations without 
any further simplifying assumption. In this approach, changes in the metric coefficients are 
directly related to changes in the tetrad null vectors. Furthermore, perturbations in the 
vectors and lead to ffist order alterations in Weyl and Ricci scalars, and also in 

the spin coefficients. The full system of linearized perturbation equations is obtained by 
expanding the complete NP equations on the plane-symmetric black hole background, and 
retaining only the ffist order terms in the perturbation functions. 

Generally, the linearized NP equations form a large system of coupled equations whose 
solutions are difficult to be found. However, for Petrov type D vacuum spacetimes, the 
problem reduces to solving a pair of decoupled equations for the gauge and tetrad invariants 
\l/o and \&4 3]. In the present case, these equations assume the simple forms 




The foregoing equations are already linearized since the Weyl scalars and \&4 vanish in 
the equilibrium state, and then they are to be considered ffist order perturbations. The 
other quantities take their unperturbed values given in Eqs. ©-(El)- 

The background symmetries, together with the perturbation equations (|7j) and (jH)), indi- 
cate that each one of the perturbation functions \l/o and \l/4 can be written as a product of 
four different functions depending upon only one of the coordinates t, r, Lp, and z. Moreover, 
due to the wavelike character of the equations, the dependence of these functions on the co- 
ordinates t, if and z is of the form exp[i(ci;)f: + k^pip + k^z)], where u is the frequency, and 

and k^ are wave numbers which may assume discrete or continuous values, depending on 




[(A + 37 - 7* + 4/i + - p) - 5*5 - 3^2]^4 = 0. 



[{D - 4p - p*)(A - 47 + p) - 55* - 3^2]^o = 0, 



(7) 
(8) 
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whether and z are compact or noncompact coordinates, respectively. The radial functions, 
in turn, are conveniently written as 

where the subscripts ±2 are related to the spin weight and to the conformal weight of the 
Weyl scalars \l/o and \l/4. 

The resulting equations for Y±2 take the following standard form: 

AV±2 + PA^Y±2 - QY±2 = 0, (9) 
where we introduced the differential operators 



dr^ ' dr? 



Here is the tortoise coordinate, defined in such a way that dr/dr^ = f{r), and 
with k"^ = kl + kl. 

The Chandrasekhar transformation theory can now be used to transform Eqs. (jH} 
into one-dimensional Schrodingerlike equations of the form 

^2^(±) ^ y{±)zi±)_ (10) 

This is done by means of the substitutions 



whose inverse transformation gets the form 



^4 ^4 



K^^^Z^^^ = -QY+2 - ^^{W^^^ + 2iu)k_Y+2, 



_^4 ^4 

-pQy-2 - 



where 

= fc^ ± 24icjM. (11) 
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In the foregoing equations, the potential functions V^"^^ are given by 



576M3 + 12k^Mr^ + k^r^ + lUM^r{k'^ + 2r' 
(Pr + 12M)2 



(12) 



= ^(A;V- 12M), (13) 

while the auxiliary functions W^"^^ are, respectively, 

^ 12M(2r3 + k'r + 4M) ^ _12M 

r2(Pr + 12M) ' ' 

The meaning of the signs (±) comes from the metric approach to the perturbation prob- 



lem, as it was performed in Ref. j37|. They are related to the parity property of the 
perturbed metric functions under the exchange —Lp (or z —z). The variable Z^"*"^ 
represents the polar (even) perturbations, while Z*^^-* represents the axial (odd) ones. 

In a certain sense, the wave equations (jlUj) can be thought of as a generalization to those 
found by Cardoso and Lemos '37|. Here the perturbations depend on both coordinates 



and z, while in the aforementioned work the analysis was restricted to ^-independent metric 
variations. This is equivalent to the vanishing of the wave number kz-, resulting in k"^ = k'^. 
However, the plane symmetry of the background spacetime renders possible to find new 
coordinates (</)', z') in such a way that one of the wave numbers, k'^ or A;^, is zero. That is 
to say, it is always possible to find a coordinate system in which the perturbations depend 
only on two space coordinates (r and or r and z). 

Another result derived from the above analysis is related to the Starobinsky constant j^|, 
|^|2 = K'^'^^K'^'^\ where i^^^^ are given by Eq. (fTTj) . The solutions of = correspond 
to algebraically special frequencies in the sense of Chandrasekhar [4^. As we shall see in 
Sec. IVIf one of such solutions, given by 

LUa = i—TTT = i T, (14) 

24M Qrl' ^ ^ 

appears in the numerical results in connection to the purely damped QNMs. The properties 
of the axial perturbations related to these algebraically special modes are analyzed in some 
detail in Sec. IVTTTl 
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IV. QNM ANALYTICAL PROPERTIES 



To solve the wave equations (fTUj) we must impose boundary conditions on the possible 
solutions. First, however, consider the case of the asymptotically flat Schwarzschild space- 
time. Gravitational perturbations lead to Schrodingerlike equations whose potentials vanish 
at inflnity and at the event horizon. Hence, close to these extremities, the solutions are 
plane waves, Z gi*'^^'*^ where the r^, coordinate ranges from — oo to +cxd. QNMs are then 
defined as solutions which are purely ingoing waves at the horizon, Z — > e"'"*'^'"* , and purely 
outgoing waves at infinity, Z —>■ e"*'^''* . These requirements are satisfied only for a discrete 
set of complex u called the quasinormal frequencies. 

In asymptotically AdS spacetimes, we keep the requirement of purely ingoing waves at 
the horizon. In the meantime, since now the tortoise coordinate takes a finite value at 
spatial infinity, the boundary condition at infinity must be changed. This can be Dirichlet, 
Neumann, or Robin boundary condition, depending on whether the field, its derivative, or 
a combination of both vanishes, respectively. Here, as in most of the literature, we adopt 
Dirichlet boundary conditions, since these lead to energy conservation in the field theory on 
AdS space (for alternative boundary conditions see Ref. jl2|). 

Having defined the QNMs, we show below that the imaginary part of the QN frequencies 
is positive for non-negative potentials. The demonstration follows the one carried out in the 
Schwarzschild- AdS case by Cardoso and Lemos |9]. Initially, we exchange the variable Z by 
= e-''^''*Z, so that Eq. ^ becomes 

where Z, and V stand for both perturbation types, the polar {Z^'^\ (f)^^\ V^^'^) and the 
axial ones {Z'^~\ (f)^~\ V^~'>). 

Near the horizon and at infinity, the solutions satisfying the QNM boundary conditions 
are given, respectively, by 



jCi/3rv)^+ 



K = ^T^7^^77^7^:7t;T^' (16) 



0OO = 5 e'"/'' sin (Vcu2 - Cs/r) , (17) 

where A and B are arbitrary constants, and Ci = V{r^)/ f{r^) and C2 = V{oo) are definite 
constants whose values depend on the perturbation type. 



Now we multiply Eq. (fT^ by 0* and integrate by parts using the asymptotic solutions 
((1^) and (fT7|). After some algebra we obtain 



/ 



dr 



*=MMr±£ (18) 

Im(co') 



Since /(r) assumes only positive values outside of the black hole, the sign of the quantity on 
the left hand side of Eq. (fT^ depends crucially on the sign of the potential V{r). For non- 
negative V, it follows that the imaginary part of the QN frequencies is positive, Q{uj) > 0. 
Moreover, it is seen from Eq. ()12j) that V^^^ is positive in the region r+ < r < cxd for any k 
and r+, and then we conclude that the plane-symmetric AdS black holes are stable against 
polar perturbations. The same is true for axial perturbations as far as the wave numbers 
are restricted by > \/3r^. On the other hand, for k < \/3r+, the potential V^~^ takes 
negative values and the theorem cannot be applied. 

Another interesting property of frequencies associated to perturbations of a plane- 
symmetric black hole comes from the particular symmetry of such a background. This 
symmetry allows us to rescale u, k and r+ by a pure scale transformation of coordinates 
t = at', (p = a{p', z = az', and r = r'/a, where a is a constant [J]. It implies that the 
frequency is a first degree homogeneous function of its parameters: u{ar^, ak) = auj{r^, k). 
In particular, we can take a = k^/k so that 

k f ko 



" VT'^^' 7 ■ ^^^^ 
Then, after computing the QN frequencies for a given /c, e. g., for k = ko, and for different 
values of r_|_, we can obtain uj{r^, k) by using the previous equation, and the values of uj for 
a fixed wave number k are naturally divided into three different regimes: large (r+ ^ k), 
intermediate (r+ ~ k) and small (r+ <^ k) black holes. In the large (small) black hole limit, 
k (r+) is negligibly small and it is expected that oc r+ {u oc k). 

We are now in a position to solve the perturbation equations (fTUj) and to find the values of 
the QN frequencies. This task, however, is more easily done through numerical procedures. 
Therefore, it is appropriate summarizing here the method used for such a purpose. 

V. THE HOROWITZ-HUBENY METHOD 



In the last decades a series of methods has been developed to numerically calculate 
QNMs of black holes. Here we briefly describe one of these methods, used in Ref. 



he 



which is of particular interest for the present work since it is suited for asymptotically AdS 
spacetimes. In order to apply such a method we rewrite here the perturbation equations in 
the appropriate form. 

By introducing the independent variable x = 1/r, Eq. (jl5|) becomes 

{x'-xl)^ + {3x' + 2tu;xl) ^ + V<P = 0, (20) 
^ dx^ ^ ax 

where x+ = l/r+ and V{x) = x\V/ fx"^, with V and / now being considered as functions 
of X. The foregoing equation has a regular singular point (x = x+) in the region of interest, 
< X < x+. Then, by Fuchs theorem, any solution of ()20|) can be written as a Frobenius 
series js^. The behavior of the solutions near the horizon may be explored by writing = 
(x — x+)°, where a is a constant to be determined. By substituting this expression into Eq. 
(PUj) and keeping terms only to the leading order, it follows that x\a{?>a + 2iujx^) = 0, which 
has the roots a = and a = ~2iuj/3r^. These roots correspond exactly to the ingoing and 
outgoing waves at the horizon, respectively jj]. Since we want only ingoing modes, we take 
a = and look for solutions of the form 

oo 

0(x) = ^a^(x-x+)'". (21) 

m=0 

Substituting Eq. (j2H) into (j2(J|) . we find expressions for the a^'s in terms of the frequency 
uj (and of the parameters k and x+). Then, imposing the boundary condition at infinity, 
0(0) = 0, we obtain 

OD 

Y^dmi^) i-^+r = 0. (22) 

m=0 

Thus, the problem of finding QN frequencies has been reduced to that of obtaining the roots 
of the infinite-order polynomial equation ()22j) . Of course, we cannot carry out numerically 
the full sum in expression (j^^ . In effect we truncate the series after a large number of terms 
and look for the zeros of this partial sum. The precision of the results is then verified by 
computing the relative variation between the zeros as we perform higher partial sums. We 
stop our search once we have a 6 decimal digit precision. For each pair of k and r_|_, we find 
a set of QN frequencies (out of the infinity of possible frequencies). These frequencies are 
labeled as usual with the principal quantum number n and ordered in a set beginning with 
the roots with the lowest imaginary parts. 
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VI. NUMERICAL RESULTS AND ANALYSIS 



We use the Horowitz-Hubeny procedure described above to compute the gravitational 
QNMs of plane-symmetric AdS black holes. As it happens to other AdS black holes, there 
are purely damped modes, whose frequencies are pure imaginary numbers, in addition to the 
ordinary modes, the frequencies of which also have nonzero real parts which are responsible 
for the oscillatory behavior of the perturbations. 

As mentioned earlier, the numerical analysis is divided into three different regimes ac- 
cording to the ratio r^/k values. That is to say, we have (i) large black holes when the ratio 
r^/k is large compared to unity; (ii) intermediate black holes if r^/k is of the order of unity; 
and (iii) small black holes when r^/k is small compared to unity. 

It is also worth mentioning here that some numerical analysis of the gravitational QNMs 



was done in Ref. [37[. However, such a study was restricted mainly to the large black 
hole regime for a specific value of k, without generalizations for a generic wave number. In 
what follows, we perform a more detailed analysis, including small and intermediate black 
holes, and extend implicitly the numerical results to all values of k. We also search the QN 
frequency values in the limit of highly damped overtones. 



A. Ordinary quasinormal modes 

To begin with, we discuss the case of regular (ordinary) QNMs for which the frequencies 
have nonvanishing real and imaginary parts, uj = uJr + ioJi. We list some numerical results of 
the lowest quasinormal frequencies for k = 2 and some selected values of r+. The computed 
quantities include both the axial and polar perturbations, and the three regimes defined 
above. The data for large, intermediate, and small black holes are shown, respectively, in 
Tables HHni and HTH 

The first feature to be noticed is that, contrary to the asymptotically flat space case, the 
axial and polar perturbations no longer have the same spectra. In fact, the isospectrality is 
only restored in the large and small black hole limits (see also Refs. IfS^). 

The displayed numerical results also allow us to obtain the k dependence of the modes 
on basis of the homogeneity property of the function a;(r+, k), as represented by Eq. (fTUIl. 
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hole regime. 





Axial 


Polar 








UOr 




1000 


1849.42 








500 


924.710 


1331 92 


924.727 


1331 92 


100 


184.948 


266.384 


184.963 


266.351 


50 


92.4820 


133.189 


92.5134 


133.124 


10 


18.5492 


26.6228 


18.6960 


26.2961 


5 


9.35491 


13.2880 


9.62439 


12.6329 


t ordinary QN frequencies for k = 


= 2 and several values of 




Axial 


Polar 


r+ 




OJi 


UJr 




2 


3.92888 


5.24698 


4.06031 


3.84425 


1.6 


3.21225 


4.15458 


3.21731 


2.78751 


1.2 


2.44666 


2.99891 


2.52608 


1.92268 


1 


2.04735 


2.21550 


2.30526 


1.55218 


0.8 


2.05189 


1.41550 


2.18361 


1.19971 


0.4 


2.09860 


0.520025 


2.10247 


0.511896 



TABLE II: The lowest ordinary QN frequencies for k = 2 and several values of r+ in the interme- 
diate size black hole regime. 



More specifically, the k ^ associated frequency can be found by means of the relation 

k 



uj{r+,k) 



2 \ k' 



(23) 



For large black holes, our results are in perfect agreement with those found in Ref. |37|. 
Axial and polar perturbations have real and imaginary parts of the fundamental frequencies 
given approximately by uJr = 1.85r_(. and uJi = 2.66r_|_. According to Eq. ()23|1 . these linear fits 
are wave-number independent. They are also equal to those from the large Schwarzschild- 
AdS black hole as should be expected, since a spherical system approaches a plane one 
as the radius goes to infinity. 
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Axial 


Polar 










1/10 


2 02727 




2 02727 




1/12 


2.02223 


0.0733694 


2.02223 


0.0733649 


1/14 


2.01866 


0.0608423 


2.01866 


0.0608402 


1/16 


2.01601 


0.0517477 


2.01601 


0.0517466 


1/18 


2.01397 


0.0448789 


2.01397 


0.0448783 


1/20 


2.01236 


0.0395123 


2.01236 


0.0395118 



TABLE III: The lowest ordinary QN frequencies for k = 2 and several values of r-|_ in the small 
black hole regime. 

For small and intermediate black holes, the quasinormal frequencies do not scale with 
the horizon radius. This is clearly shown in Fig. where we plot the lowest axial and polar 
frequencies as functions of r_|_, for k = 2 and 0.1 < r_|_ < 4 . The exact way the curve uj{r^) 
deviates from the linear behavior depends on the perturbation parity. Within a certain 
accuracy, the fundamental frequency starts deviating from the line uj = (1.85 + z2.66)r+ 
around r+ = k/2 for axial modes, and around r+ = 5k/ 2 for polar modes. As the ratio 
r+/fc goes to zero, we have ujr ^ k and oji 0, independently of the perturbation type. In 
fact, the dependence of u upon k in the limit of small black holes may be obtained using 
Eq. fll9|) which gives uj{r^, k) ^ uj (0, k^) k/k^ = const, x k . The numerical analysis shows 
that for the fundamental mode the proportionality constant is equal to unity. The relation 
of this result to the pure AdS modes is investigated in Appendix A. 



B. Purely damped modes 

Axial perturbations of AdS black holes are known to have special quasinormal modes 
which are purely damped. In the present case, the values of the associated frequencies, 
UJ = ius, are listed in Table HVl for k = 2 and some selected values of r+. Once again, the 
homogeneity property of uj{r^, k) allows us to obtain the k dependence of the frequency on 
the basis of relation (1^ . 

n 

For large black holes, our numerical results agree with those presented in Ref. fST]. To 
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AXIAL 




FIG. 1: The frequency of ordinary QNMs for axial and polar perturbations, cu = uJr + iooi, as 
a function of the horizon radius r-|_ for small and intermediate black holes. The dots represent 
numerical results, and the solid lines are linear fits. On the left are the real parts (w^) with the 
linear fit 1.85r+, on the right, the imaginary parts (wj), with the fit line 2.66r4.. 

a remarkable accuracy, the k = 2 pure imaginary mode is well described by the formula 
LOs = 4/3r_|.. The extension of this fit to any k leads to uJs = k'^/3r^. A similar result was 



already obtained for the Schwarzschild-AdS spacetime by Cardoso et al [31[. They found 
ujg = {I — !)(/ + 2)/3r+, where / is the angular momentum of the perturbation. The above 
two formulas for uo^ are just fits to the numerical data, but we were able to find an exact 
analytical proof for both of them by expanding ujs in a power series of l/r_|_. This study is 
the subject of the next section. 

For intermediate size black holes, the frequency uOg increases faster than /i;^/3r+ as the 
horizon radius decreases (see Fig. |2I). Another special feature of this regime is related to the 
behavior of the axial potential V^(~) as a function of the coordinate r. As we go from larger 
to smaller horizon radii in the plane r+ x fc, the straight line k = v^'"+ is the boundary on 
which the axial potential (fT!?|) becomes a non- negative function in the region r+ < r < oo. 
In other words, V'^^\r) > at all points outside of a black hole perturbed by an axial 
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4/3r+ 


r+ 




r+ 




8/3rl 


100 


0.0133342 


0.0133333 


2.0 


0.733037 


0.9 


3.66330 


3.65798 


75 


0.0177794 


0.0177778 


1.8 


0.835965 


0.7 


7.77447 


7.77454 


50 


0.0266704 


0.0266667 


1.6 


0.978431 


0.5 


21.3333 


21.3333 


25 


0.0533639 


0.0533333 


1.4 


1.19459 


0.3 


98.7654 


98.7654 


10 


0.133786 


0.133333 


1.2 


1.58723 


0.2 


333.333 


333.333 


5 


0.270361 


0.266667 


1.0 


2.64636 


0.1 


2666.67 


2666.67 



TABLE IV: The pure imaginary QN frequency for k = 2 and some selected values of r^. The 
formula 4/3r-|. is a fit valid for large black holes, while 8/3r^ is the magnitude of the algebraically 
special frequency tOa- 

8 
6 

m 

4 
2 



0.8 1 1.2 1.4 1.6 1.8 2 

r+ 

FIG. 2: The frequency cOg of the k = 2 purely damped mode for intermediate and small black 
holes as a function of the horizon radius (the dots represent numerical results). The dashed line 
is u!s = 1/3?"+, and the solid line is the algebraically special frequency |a;a| = 8/3ri^. At point 
P = (2/ \/3, Vs), the algebraically special frequency line crosses the purely damped QN frequency. 

mode with k > ^/Sr^^. A consequence of this for the numerical calculations is that, at 
= k/ \/3, our root-searching program based on the Horowitz-Hubeny method fails to find 
a pure imaginary QN frequency. The problem occurs when r_(_ approaches the point where 
the function u = iujs{r^) would intercept the algebraically special frequency Ua = ik'^/6r^. 
This is indicated by the point P in Fig. |2l where we plot Ug as a function of r+ for k = 2. 
The numerical data show that tUg > \uJa\ for all r_|_ > k/^/S, and Ug < \uja\ for all r+ < k/y/S. 
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n 


Axial 


Polar 










1 


1 RA OAS 




1 84 963 




2 




491.641 


31 fi 1 59 


491 583 


3 


446.438 


716.753 


446.482 


716.671 


4 


576.528 


941.807 


576.586 


941.700 


5 


706.536 


1166.84 


706.608 


1166.71 


20 


2655.34 


4541.91 


2655.62 


4541.40 


40 


5253.44 


9041.92 


5254.09 


9041.02 



TABLE V: QNMs corresponding to = 2 axial and polar perturbations of a large plane-symmetric 
AdS black hole, with r+ = 100. 



For small black holes, the purely damped frequency iug approaches the algebraically 
special frequency uja from below. In fact, the difference between these frequencies is smaller 
than 10~^ for r^/A; = 0.3 and decreases as r^/A; — > 0. This conclusion has again its spherical 



counterpart, as it can be seen from the numerical results found in Ref. 



3l| . More details 



on the algebraically special frequency and on its non-QN character are given in Sec. IVIIII 



C. Highly damped modes 

For the moment, we present only the computed frequency values of the lowest QNMs. 
They are the most important modes for the AdS/CFT correspondence, since these give 
the thermalization time scale on the boundary field theory. However, as soon as one is 
interested in Hod's conjecture js^, which relates the classical QN frequencies with the 
quantum spectrum of black holes, the most relevant modes are the highly damped QNMs. 
Indeed, according to such a conjecture, the real part of the QN frequencies, in the limit 
of high overtones, gives information about the quantum of the black hole area. In the 
following, we present the QN frequency values for modes in the regime of high principal 
quantum numbers n. Since we can distinguish the purely damped modes as belonging to 
a special family, we simply label them with n = 0, and begin to label the fundamental 
ordinary modes with n = 1. 
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n 


Axial 


Polar 










1 






9 30596 




2 




4 74580 


3 961 37 


3 791 63 


3 


4.65907 


7.05251 


4.39306 


5.96477 


4 


5.95653 


9.32552 


5.59513 


8.21840 


5 


7.24999 


11.5884 


6.83055 


10.4733 


20 


26.6726 


45.3810 


26.0793 


44.2560 


40 


52.6238 


90.3921 


44.2193 


75.7648 



TABLE VI: QNMs corresponding to = 2 axial and polar perturbations of an intermediate size 
plane-symmetric AdS black hole, with r+ = 1. 



n 


Axial 


Polar 


OJr 








1 


2.05667 


0.214983 


2.05678 


0.214687 


2 


2.17823 


0.576255 


2.17845 


0.575676 


3 


2.32373 


0.974834 


2.32406 


0.974002 


4 


2.48740 


1.39433 


2.48783 


1.39325 


5 


2.66559 


1.82658 


2.66613 


1.82528 


20 


6.04295 


8.60462 


6.04536 


8.60001 


40 


11.0316 


17.6657 


11.0366 


17.6567 



TABLE VII: QNMs corresponding to A; = 2 axial and polar perturbations of a small plane- 
symmetric AdS black hole, with r_|_ = 0.2. 

To study the QNMs asymptotic behavior, we computed the first 40 ordinary QN frequen- 
cies associated to = 2 axial and polar perturbations of plane-symmetric AdS black holes. 
A few representative results for some modes of large, intermediate, and small event horizons 
are displayed in Tables El IVIl and IVIIl respectively. 

On the basis of this numerical data, we obtained linear fits by taking overtones from 
n = 21 to = 40 and employing the discrete least squares approximation method. Together 
with the coefficients of the linear approximating functions, we calculated the error involved 
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Axial 


Polar 




a 




a 




100 


57.2447 + 41.9027i 


129.905 + 225.000i 


57.0927 + 41.7360i 


129.925 + 224.982i 


1 


0.717273 + 0.371668i 


1.29763 + 2.25053i 


-7.58881 - 14.2636i 


1.29504 + 2.25074i 


0.2 


1.02468 - 0.451644i 


0.249900 + 0.452989i 


1.02451 - 0.451813i 


0.250029 + 0.452767i 



TABLE VIII: Linear fits for the frequencies of high overtones of axial and polar perturbations for 
black holes with r_|- = 100, r+ = 1, and r+ = 0.2. 

in each fit process and established 1% as the upper bound error. This means that, in the 
regime of high overtones {n oo), the QN frequencies are evenly spaced with real and 
imaginary parts given approximately by linear functions as described below. 

The linear fits depend upon the size of the horizon radius, and are given by functions of 
the form 

ujn = a + Pn , (24) 

where the coefficients a and P refer to both the real and the imaginary parts of the frequency, 
and are given in table IVmi Their dependence upon the black hole parameter r_|_ and on the 
wave number k is governed by Eq. ()19|1 . 

As should be expected, the results found for large black holes reproduce those obtained 



in Ref. [31[ for Schwarzschild-AdS black holes with a large event horizon. A special feature 
of the regime is that axial and polar perturbations have essentially the same QNM spectra, 
which appear in the first line of Table IVIIII In the other regimes, however, the linear 
functions are very different from those of the spherical case. 

The angular coefficients of the linear functions seem to present a kind of universal be- 
havior. In fact, as seen from the P coefficients in Table IVIIH the results for both of the 
perturbation types, and for each horizon radius, are very close to each other. The spacings 
between the frequencies per unit horizon radius of consecutive modes are essentially inde- 
pendent of the black hole size. This is seen more clearly in Table HXl where the values of the 
ratio {uJn+i — ^n)/'r+i calculated for each regime, are shown. The data confirm that such a 
ratio, for the calculated wave number = 2, is approximately constant. Additionally, using 
this result together with Eq. ()23|) it may be shown that, in the high overtones regime, the 
frequency per unit horizon radius spacing is not only independent of the black hole size and 
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(^n+l ~ 




r+ 


Axial 


Polar 


100 


1.29905 + 2.25000i 


1.29925 + 2.24982i 


1 


1.29763 + 2.25053i 


1.29504 + 2.25074i 


0.2 


1.24950 + 2.26495i 


1.25015 + 2.26384i 



TABLE IX: The fractional spacing of frequencies among consecutive high overtones of axial and 
polar perturbations of black holes with r_|_ = 100, r+ = 1, and r+ = 0.2. 

of the perturbation parity, but it is also independent of the wave-number value. The real 
parts for small black holes are about 4% lower than for intermediate and large black holes. 
We notice that this difference decreases as we search for overtones with principal quantum 
numbers much larger than n = 40, but we have not investigated this more precisely because 
the numerical calculations become very slow for overtones with n significantly larger than 
40. 

VII. THE LARGE BLACK HOLE LIMIT OF THE PURELY DAMPED QNM 

Although numerical procedures need to be used to compute the QN frequencies, it is 
possible to study analytically the special case of axial perturbations of very large black 
holes. Here we show in some detail the process of finding an asymptotic expression for the 
frequency of the purely damped modes. The analysis is performed in such a way that it 
holds for both the plane-symmetric and the Schwarzschild-AdS spacetimes. 

Our starting point is the Horowitz-Hubeny approach described in Sec. |3 There, the QN 
frequencies of axial perturbations are associated to solutions of the equation 

oo 

^a„M (-x+r = 0. (25) 

m=0 

Hence, we need to obtain an analytical expression for the sum of this series in order to find 
exact frequency values. Such a task is very difficult to be carried out for a general but 
it becomes feasible for the large black hole limit of purely damped modes. 

The coefficients 0^(1^) appearing in Eq. (^3)) are the Frobenius coefficients of the power 
series solution to either Eq. ()2()|1 written for the axial functions, or the corresponding 
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equation in the Schwarzschild-AdS case (see Ref. They are given in terms of the 

frequency, and of the other physical parameters of the problem, by the following recurrence 
formulas: 



ai(-x+)i = AFiao(-x+)°, (26) 

a^{-x+r = — [Amam-i{-x+r-' - i?„a„„2(-x+)"-'] . (27) 
m 

The last equation is valid for all m > 2 and we have shortened the equations by introducing 
the quantities 

Am = 3(m2 - m - 1) + [q{q + k) + {2m^ - 2m - 3)k] x^, (28) 

Bm = (m^ - 2m - 3) [1 + kxI] , (29) 

Fm = [m(3 + Kxl) + 2iujx+] , (30) 

where q = k,l and k = 0, 1 depending on whether it is considered a plane or spherical hole, 
respectively, and now the foregoing expressions hold for any m > 1. 

In the large black hole regime, the wave number (or the angular momentum) of the 
perturbation is negligible and the frequency scales as ~ l/a;+. In terms of the variable 
x+, it implies that u has a first order pole at x+ = 0. Therefore, uj{x^) has a Laurent 
series about x+ = 0, which is valid over the range < x+ < -R, for some R 0|. Such an 
expansion can be conveniently written as —2ix^u! = Ylm=o ^mX^, where are constants to 
be determined. However, not all of the elements 6^ are independent, since part of them can 
be eliminated on the basis of the parity property of the function uj{x+). Independently of 
what results from the sum on the left-hand side of Eq. (j^ . the form of the quantities A^, 
Brni and ensures that the product x^uj is a function of x^. This means that uj changes 
its sign under the transformation being an odd function of Xj-.. Accordingly, the 

series representing x+uj should be an even function of x+, which means that only the even 
powers of x+ should survive, and then one should have 62n+i = for all > 0. This result 
allows us to write 

oo 

-2ix+uj = J2 Cm^", (31) 

where z = x'^ and Cm are constants to be determined. 

The series representation of 2ix^uj opens the possibility of Taylor expanding the terms 
<^m{—x+)"^ about z = 0. Before doing this, we notice that the value of Oq is arbitrary, since 
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all the terms in the series of Eq. ()25|) are proportional to this element and the full sum is 
equal to zero. For simplicity, we choose here oq = 1. Then we substitute Eq. into 
Eq. (jHU)) and expand (which depends on x+ and also on uj) about z = When this 
expansion is substituted back into Eqs. (j26|) and (j27j) . we obtain a set of Taylor series for 
the elements am(— x^)'" (one of such series for each m), which may be represented by 



j=0 



where the elements rj^ are functions oi q, k, Cq, Ci, Cj. Substituting the previous expan- 
sions into Eq. (^3)) and collecting terms of the same power in z, we find 

[1 + r7°(co) + V2ico) + •••] + Hico, ci) + vlico, d) + ...] z + ... = 0. (32) 



Yet the left-hand side of Eq. ()32|) is a power series in the variable z. It implies that, 
in order for the equation being satisfied, all bracketed terms should vanish. Therefore, 
the zeroth order coefficient provides an equation for cq, the first order coefficient, together 
with the earlier calculated value of Cq, provides an equation for ci, and so on. Here we 
begin with the computations that lead to Cq. The first two relevant terms for this task 
are ?7° = — 3/(3 — cq) and = — 3co/[2(3 — co)(6 — cq)]. Moreover, it can be seen from 
Eq. fl29|) that is identically zero, and as a consequence the term 03 (— x_|_)^ differs from 
a2 (— x+)^ only by a multiplicative factor. This implies, in particular, that the element r)^ 
is related to by r]^ = 5/(9 — co)?72. Such a relation, together with the formulas (j26|l and 
(I27p . implies that all the subsequent elements are also proportional to r/g- More specifically, 
the mth coefficient (for m > 2) can be written as r]^ = 7^(00)7725 where 7m(co) is given by 
the recurrence formula 

7m(co) = —TTT— T [3(m^ - 171 - l)7„^_i(co) - (m^ - 2m - 3)7^-2(00)] , 

m[6m — Co) 

with 72(co) = 1 and 73(00) = 5/9. 

Taking into account the last results, the ^-independent term of Eq. (fH^ gives rise to the 

following equation for Cq: 

1 3 



Co 



(3 -Co) 2(3-co)(6-Co) 



m=2 



0. (33) 



Clearly the above equation has the special solution co = 0, in addition to the roots that 
arise as zeros of the square-bracketed term. The existence of this special solution explains 
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why the purely damped modes do not scale with the horizon size, as it happens to ordinary 
modes in the large black hole limit. In fact, the numerical computations show that purely 
damped modes are the only QNMs which in first approximation are not proportional to r_|_. 
This fact suggests they should be associated to the root Cq = 0, and then it follows from Eq. 
(|31|) that for large horizon radii the corresponding frequencies go as l/r+. On the other hand, 
the numerical results also indicate that the ordinary axial QNMs are proportional to r+ (in 
first approximation). From this one concludes that the remaining roots that are obtained 
from the zeros of the term among brackets in Eq. PHj) correspond to the asymptotic values 
of the function 2ix^ uo for these ordinary QNMs. 

It is seen from the above analysis that, in order to calculate the leading term of the large 
horizon limit for the purely damped mode frequency, we must take Cq = and calculate the 
next parameter ci. In such a case, the coefficients vi\ and rj\ reduce to vi\ = [q{q + — 
2k, — Ci]/3 and r]l = — Ci/12, while the other elements (for m > 3) assume the form 
vln = 7m (0)772- Substituting these results into the first order coefficient of (l32|) we obtain an 
equation that can be solved for ci, yielding 

_ 4 [q{q + 2k] 

where we have defined 

CO 

m=0 

The problem now is summing the series that leads to 7. Fortunately, we succeeded in 
summing up this series (see Appendix B) and the exact result is 7 = 2. With this, the 
coefficient Ci becomes 



2 

ci = ■^{q-q-){q-q+), 



where g_ and q+ are given by 



-K±y/K'^ + 8k . . 

Q± = 2 ■ ( ^ 



Once we have obtained an expression for ci, the asymptotic formula of the pure imaginary 
QN frequency follows directly from Eq. (jHT|) . In fact, considering a plane-symmetric AdS 
black hole, for which q = k and k = 0, we have q^ = q^ = 0, and the interesting result 
follows: 

u = i- 



3r_| 
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Alternatively, for a Schwarzschild-AdS black hole one has q = I and k = 1, which implies in 
g_ = 1 and g+ = —2, leading to 

. (^-l)(^ + 2) 
u = I . 

Both of the last two formulas for the purely damped frequencies are exact results at the 
asymptotic limit of large black holes, in perfect agreement with the approximate fits found 
through the numerical calculations. 



VIII. ALGEBRAICALLY SPECIAL MODES 



The numerical results of the QN frequencies revealed a close relation between the al- 
gebraically special modes and the purely damped QNMs. For a fixed wave number, the 
imaginary frequency uj = iUg approaches the values of the special frequency Ua as the hori- 
zon radius goes to zero (see the comments at the end of section UTT]) . There also seems to 
exist a meeting point P between the two frequencies for a special combination of parame- 
ters such that y*-^''(r+)//(r+) = 0, as it can be seen from the numerical results shown in 
Fig. 121 Furthermore, the relation between the purely damped modes and the algebraically 
special frequency in the small horizon limit is not restricted to the plane-symmetric black 
hole case, but it also appears in axial perturbations in the Schwarzschild-AdS spacetime. 
All these facts render it convenient to analyze here the QN character of the algebraically 
special perturbation modes. The frequency function of these modes for the plane-symmetric 
black hole is given by Eq. (|T^ . while the corresponding expression in the spherical case is 
found in Ref. In a general form, the frequency Ua can be written as 

'^'^ " ' 6(1 + Kxl) 

where again q = k,l and k = 0, 1, as it is considered a plane or spherical black hole, 
respectively. 

To verify if the algebraically special modes satisfy the QNM boundary conditions we need 
to obtain the solution (pi \x) associated with the algebraically special frequency Ua- Since 
this function is a particular case of the general power series solution ()21|) . we may use Eqs. 
(j2ni) and (j27j) to determine the coefficients for u = Ua- Once again we may take = 1, 
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so that the coefficients oi and 02 are given, respectively, by 

3(1 + KX^) — q{q + k)x1 



x+ (3 + Kx^ + 2iijx+) 
x\q{q + n)[q{q + n) — 2k] + 6iuj{l + kx^) 



Ax-^{3 + Kx^ + iu;x_|_)(3 + nx^ + 2za;a;+) 



It is seen from the last equation that 02 vanishes for the frequency uj = uja- Moreover, since 
all the other elements (for m > 3) differ from 02 only by a multiplicative factor, we have 
03 = 04 = ... = 0. The solution for (pi \x) then reduces to 



We know from the Horowitz- Hubeny method that the function ()36|) represents an ingoing 
wave at the event horizon. Here the problem is verifying if (pi \x) satisfies the Dirichlet 
boundary condition at the spatial infinity. Imposing the condition (pir\o) = 0, we find an 
equation that is satisfied only if g = g.,. or g = where g+ and g„ are given by Eq. (j35j) . For 
a plane-symmetric AdS black hole, this means we need a vanishing wave number in order to 
have a QN frequency equal to Ua- In this case, however, the resulting perturbation functions 
would be independent of time, which could not represent gravitational waves. Similarly, 
for a Schwarzschild-AdS spacetime we need to have / = —2 or / = 1. The first root is a 
physically forbidden value, while the second one leads to a static perturbation, which only 
adds a small rotation to the black hole js^. On the basis of these results, we conclude that 
uJa cannot be a QN frequency, because the associated wave function (pi ^ (x) does not satisfy 
the Dirichlet boundary condition at infinity. 

Although the numerical calculations show that in the small horizon limit u = iuj^ is 
very close to Ua, the above result proves that the frequency function associated with purely 
damped modes is not exactly equal to the algebraically special frequency. The solid line in 
Fig. 121 does not coincide with the dots, which represent QNMs, not even in the small black 
hole regime. There is, however, a special point P where uJa seems to coincide with iug. The 
answer to the question whether Ua corresponds to a QN frequency at that point or not is 
given in the following. 

In spite of the generality of the analysis performed above, it fails to test the QN character 
of UJa in the case where V^^^''(r_|_)//(r+) = 0. This situation happens for a special family of 
parameters satisfying the condition 




(36) 



q{q + K) 




(37) 



24 



In this case the roots of the Frobenius index equation, which are a = and a = 2ia;//'(r_|_) = 
1 (/' = df/dr), differ by an integer number and the Horowitz-Hubeny method is not vahd. 
We can circumvent this problem by writing exphcitly the axial version of Eq. (|2(Jj) for u = Ua 
and parameters satisfying Eq. (jSTj). By doing this in such a way that the analysis holds for 
both the plane and the spherical AdS black holes, it is found that 



[(1 + + x^x + x^] 



d wa 



dx"^ 



+ [3(1 + KxX)x + (3 + Kx\)x+] 



dx 



-3(1+Kxi 



= 0. 

(38) 

Here x = x+ is a regular point of the above differential equation, and the solution for (pi ^ 
can be expressed Taylor series. 



[X 



Y^dr, 

r?i=0 



( X X -\ 



Substituting this expression into Eq. (jHHj) . we find as usual that all coefficients dm (for 
m > 2) can be expressed in terms of the first two, d^ and di. The general solution of Eq. 
(I38|) then assumes the form 

3(l + /txl)^ /?„ 



0i \x) = do 



2(3 + Kx^ 



+ di 



+ m=2 
2 \ oo 



-x+y 



[X — Xa 



^ ^ (3 + 2kx_^_) \^ (3m ^ 
[X - x+) + , ._2^ (-0;+)™-!^^ ~ 



(3 + KX^ 



(39) 



where the parameters jSm (for all m > 4) are given by the recurrence formula 

(m — 3)(m + 1)(1 + 



_ (3 + 2nxl) 

Pm /o I 2 \ Pm—l 



-A 



'm-2, 



(3 + Kx^) ' m(m - 1)(3 + ra^) 
with (32 = 1 and /^g = (3 + 2/€x^)/(3 + kxI). 

The first term among brackets in Eq. ()39|) represents a perturbation function which is 
constant at the event horizon and corresponds to an ingoing wave, while the second term 
represents an outgoing wave, since it goes as (x — x+) in the limit x — ^ x+. To satisfy the 
QNM boundary condition at the horizon, we must take di = and keep only the first term 
of the general solution (jHEI)- In addition, without loss of generality we may choose do = 1. 
Finally, we must verify if the function 0i~^(x) satisfies the Dirichlet boundary condition at 
the spatial infinity x = 0. By substituting do = 1, di = 0, and x = into Eq. it follows 
that 

3fl + Kx^ ) °° 



2(3 + Kx? 



m=2 
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0oo(l, X+) 



8 10 



FIG. 3: The values of (j)ooi^, function of the angular momentum /. The dots are numerical 

results, while the dashed line is ^oo(0,x+) ~ 3.29. 

For plane-symmetric AdS black holes, our numerical calculations show that (/)oo(0,a;+) ~ 
3.29, independently of the horizon radius. On the other hand, we obtain different values for 
0oo(l; x^) as we change the Schwarzschild black hole size. Another way to see this variation 
is in terms of the angular momentum, since x+ and / are coupled through Eq. (jSZ|)- In Fig- 
we plot the numerical values of 0oo(O, x^) (dashed line) and of (pooi^, x^) for / > 2 (dots). 
As the angular momentum goes to infinity, which means x+ 0, the results found for the 
spherical case tend to the value of the plane black hole. Since (j)oo{i^,x+) is always different 
from zero, we conclude that the algebraically special modes are not QNMs, not even for the 
special parameters that lead to V^~\r^) / f{r^) = 0. For the plane-symmetric black hole 
and for k = 2, this corresponds to the point P in Fig. |21 where the lines representing the 
purely damped QN and the algebraically special frequencies tend to intercept each other. 
Among other things, this result implies that the QNM spectrum of some black holes does 
not present a particular pure imaginary frequency, for those specific modes whose wave 
numbers (or angular momentum) satisfy Eq. (jHTj) . For instance, an axial quadrupole (/ = 2) 
perturbation is not able to excite a purely damped QNM in a Schwarzschild-AdS black 
hole with r+ = 1. Whether this conclusion has some physical implication to the AdS/CFT 
conjecture or not is still an open question. 
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IX. FINAL COMMENTS AND CONCLUSION 



We have studied the gravitational quasinormal modes of plane-symmetric anti-de Sitter 
black holes. These modes are responsible for the late time behavior of small gravitational 
disturbances. The wave equations governing the axial and polar perturbations were derived 
by a procedure based on the Newman-Penrose formalism and on the Chandrasekhar trans- 
formation theory. This approach allowed us to treat equally the two spacclike coordinates ip 
and as well as provided the algebraically special frequencies associated to the gravitational 
waves radiated by the black hole. 

Using the Horowitz-Hubeny method, we have computed the quasinormal frequencies for 
k — 2 and several values of r+ in the large, intermediate, and small black hole regimes. 
The wave-number dependence of the modes was obtained on the basis of the homogeneity 
property of the function uj{r^,k). For large black holes, the ordinary frequencies scale 
linearly with r_|_, and the pure imaginary frequency goes as l/r_|_. For small black holes, the 
fundamental modes present a;^ — > A; and cui ^ 0, while uj — iujs approaches the algebraically 
special frequency uja — ik^/Qr^. In the limit of highly damped overtones, we found the QN 
frequencies are evenly spaced, with the spacing per unit of horizon radius being a constant 
which is independent of the wave-number value, perturbation parity, and black hole size. 

In relation to stability, we have shown analytically that for a non-ncgativc potential the 
imaginary part of the frequency is positive definite for every mode. From the numerical 
results it also follows that all QN frequencies have a positive imaginary part, including 
the cases where the potential takes negative values. These results conjointly imply that 
plane-symmetric AdS black holes are stable against gravitational perturbations. 

To study the large black hole limit we have developed an analjd^ical method based on 
the expansion of the frequency in powers of l/r_|_. Application of this method to axial 
perturbations showed that there is a special frequency for which the ratio uj/tj^ goes to 
zero as r+ — > oo. The first nonvanishing term is ujs — ik'^/3r+ for the plane hole and 
LOs — i{l — + 2)/3r+ for the Schwarzschild hole. It is important to observe that these 
purely damped modes are particularly long living for large black holes, which is exactly 
the relevant regime for the AdS/CFT correspondence. If the above asymptotic formulas 
could be derived on the CFT side, it could provide an important quantitative test to the 
Maldacena's conjecture. 
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As a last remark, it should be noted that the discussion of the physical interpretation 
of the modes with different wave-number values remains. Preliminary results show that 
the zero wave-number axial perturbations yield only small rotations on the system, while 
;he polar ones lead to an increase in the mass and the production of gravitational waves 
5^. This issue, however, must be investigated in more detail and shall be considered in our 
future work. 

Acknowledgments 

We are indebted to V. Cardoso for suggestions and for a critical reading of the manuscript. 
We also thank J. P. S. Lemos for stimulating conversations. One of us (A. S. M.) thanks the 
Conselho Nacional de Desenvolvimento Cientifico e Tecnologico - CNPq, Brazil, for a grant. 

APPENDIX A: THE PURE AdS PLANE- SYMMETRIC GRAVITATIONAL 
MODES 

In this appendix we discuss briefly how to find the modes of gravitational perturbations 
in pure plane-symmetric AdS spacetimes. This is to be compared to the small black hole 
limit QNMs, as calculated numerically in Sec. IVII 

The relevant wave equations for axial and polar perturbations may be obtained from Eq. 
fllO|) by assuming M = wherever it appears. For instance, the potentials V^"^^ reduce to 
the same constant value 



(Al) 



The resulting equations are then 




(A2) 



whose solutions in terms of the coordinate r 



l/r=K are of the form 




2 ; 2 

UJ = k , 




(A3) 



where a^^^ fe^^-*, c*^^-', and d^"^^ are integration constants. 
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The solution functions for cu^ 7^ k'^ are analytic in the whole range of the radial coordinate 
which is of interest for comparison to the small black hole QNMs, < r < 00. Hence, the 
only condition uj must satisfy is 

00^ > (A4) 

in order to avoid the wave growing exponentially in the limit r ^ 0. 

However, additional restrictions to these modes arise by imposing the usual boundary 
conditions for QNMs in AdS spacetimes: no wave can come out from the black hole horizon, 
which is located at r = r+ — > 0; and no wave can propagate to spatial infinity, r 00. By 
substituting these conditions into Eqs. ()A3|) it follows that a^^-* = and 6'-^^ = 0, which 
means waves whose frequencies satisfy the relation uj^ 7^ k'^ do vanish everywhere. 

We are then left with the solutions for which u;^ = k"^ only. In this case there is in 
fact no wave in the radial direction, so that the solutions = c^^^ /r + S'^\ for r > 0, 
represent waves traveling along the spatial directions orthogonal to the r direction. In other 
words, there are no waves coming out from the horizon nor propagating to spacelike infinity. 
Considering these as the plane-symmetric pure AdS modes, their frequencies, given by 

a; = /c , 

agree with the result found numerically for small black holes. 

APPENDIX B: SUMMATION OF THE SERIES FOR MODES OF LARGE 
BLACK HOLES 

Here we show how the series presented in Eq. (I34j) can be summed up exactly. Since it 
is an alternating series, it is important not only to prove its convergence, but also to test 
it for absolute convergence. By mathematical induction, we find that the series can be 

written as a combination of seven terms 

00 00 00 00 00 00 

7 = 2+^a^+^a^+^a^+5^<+ + (Bl) 

m=l rra=l m,=l m=l m=l m=l 
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where the elements = 1, 2, 6) are given by 



33™+2(6m + 5)' ' ' 33'»(6m + l)' 

"m = (-1) o^^ir... , , «m = (-1) 



33™(6m + l)(3m+ 1)' ^ ' 33™+2(3m + l)(2m + 1) ' 

^5 = 6^ + 5 ^6 ^ , +11 

- ^ > 33-+2(3m + 2)(2m+ 1)' ^ ' 33'»+2(3m + 2)(6m + 5) ' 

Then we are left with a sum of six series, and by proving the absolute convergence of each 

one of them, the absolute convergence of the whole series ^m7m(0) is also proved. In order 

to do that, we calculate the limit 



lim 



m^QO \a\^\ 3^ 

and verify that it is smaller than unity for any z, and then by applying the D'Alembert ratio 
test we conclude that every individual series is absolutely convergent. This result implies 
that the final sum for 7 is independent of the element ordering in the alternating series. In 
particular, we can add the terms over all for i = 1, 2, 6, before performing the sum 
over m. By doing that, we find the simple result = for any value of m, and hence 

the sum over m gives zero. Therefore, the sum of the series (jBljl is 7 = 2. 
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